A simple model for the kinetics of packaging of DNA in to a capsid against an 

external force 
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We propose a simple model for the kinetics of packaging of viral DNA in to a capsid against an 
external force trying to prevent it. The model leads to a Butler- Volmer type dependence of the rate 
of packaging on the pulling force F. 



Recently, in a very interesting experiment]^, the rate 
of packaging of a long viral DNA in to a capsid by means 
of the portal complex, under the influence of an external 
force trying to prevent the packaging has been studied 
(see Fig. |l|) . The force- velocity {F — v) curves show that 
packaging rate decreases even for small forces, implying 
that the rate determining step in the packaging is af- 
fected by the externally applied force. Using a Kramer's 
type model, involving steps that involve thermal acti- 
vation over a barrier, the experimental data was fitted 
assuming two force dependent steps. The first involves 
a conformational change of only 0.1 nm and is the rate 
determining step for small forces. At higher forces, the 
velocity is found to decrease more sharply and it was 
suggested that the second step, which has a larger con- 
formational change associated with it, becomes the rate- 
limiting step. In the following, we present a simple the- 
oretical model for the F — v curves, which seems to fit 
the data reasonably well. In our model, there is only one 
step in the activation process, and we include the possi- 
bility of the reverse process too. Interestingly the model 
leads to a formula analogous to the Butler- Volmer equa- 
tion of electrochemistry. A mechanism consistent with 
the structure of the packaging motor, determined from 
X-ray crystallography has been suggested Accord- 
ing to this, the packaging involves successive firing of five 
ATPases in the portal complex. The hydrolysis of one 
ATP molecule leads to a movement of two base pairs of 




the DNA in to the capsid. To analyze the process, us- 
ing the methods of physical kinetics, we need to have an 
idea of of the potential energy surface for the processes, 
and of the transition state. Here, it is most convenient 
to take the reaction co-ordinate to be the length of the 
DNA molecule inside the capsid, which we shall denote 
as X. Every increase in x by a{= 6. 8 A) involves the hy- 
drolysis of one ATP molecule. This hydrolysis is coupled 
mechanically to the movement of the DNA, through the 
ATPase. Therefore, as x changes by a, there is a lowering 
of net free energy of the system. 
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FIG. 1: Packaging of the DNA in to the capsid against an 
external force F trying to prevent it 



FIG. 2: Free energy profile for the packaging, x denotes the 
length of the DNA that has been packaged. The packaging of 
a length a corresponds to going from A to B (or B to C) and 
involves the hydrolysis of one ATP molecule. As more and 
more DNA gets packed in the oscillations in the free energy 
profile would decrease and the profile would get flatter and 
flatter. 



A free energy profile for the process, against the reac- 
tion co-ordinate ix) is shown in the Fig ^. Each minimum 
(A,B,C) in the curve is an initial state for the packaging 
of two more base pairs, and has one ATP molecule, ready 
to be hydrolyzed. In going from A to B, the ATP is hy- 
drolyzed to ADP (leading to a lowering of free energy) 
and two base pairs have been packed in. At B, the ADP is 
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got rid off, and another ATP is attached (not to the same 
ATPase, see ref. [D) and then the system goes over from 
B to C. The experiments are done under saturating con- 
centration of ATP. Hence ATP attachment is not the rate 
determining step. The transition state is at a distance 
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FIG. 3: Packaging rate v against the force. The best fit 
was obtained using the non-linear fit package of MATHE- 
MATICA and gave v^f = 103.416p/s, v^t, = 0.0129&p/s and 
Aa; = 0.82lA. The data were taken from the figure of 3b 
of reference Q and corresponds to a packing of 1/ 3'^'* of the 
capsid 

of Aa; from A and therefore, in presence of an external 
force F, to reach it one has to do an extra work of F l\x. 
Therefore the rate of the forward process, in presence of 
the force is given by fo/ exp(— i^Ax/fcT), where I'o/ is 



the rate in the absence of the external force. The rate of 
the backward reaction is, in a similar fashion, given by 
i/obexp(i^(a — Aa;)//cT), where a = 6.8A. Hence the net 
rate is ^ = fg/ exp(— Faa/ZcT) — i/q;, exp(Fa(l — a)/kT), 
where a = Ax /a. Using this expression, we have fitted 
the F — V curve in ref. and the result is shown in 
Fig. ^. It seems to fit the data fairly well. Ref. [Q. 
assumes two different rate determining steps, operating 
at different force strengths. Further, they have four un- 
known parameters which are used for fitting. In com- 
parison, in our mechanism, there is only one (reversible 
step) and the sudden decrease in the velocity, observed 
around 45 pN is the result of the rate of reverse step be- 
coming comparable to that of the forward process. It has 
only three unknown parameters. Interestingly, the above 
equation is very similar to the Butler- Volmer equation 
of electrochemistry!^, in which the force F plays a role 
similar to the overpotential, and a is the analogue of the 
symmetry factor. The equation can be used to calculate 
the force at which the packaging velocity would become 
zero and we find it to be 55.4 pN. 
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